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Abstract 

The general properties of two-dimensional generalized Bessel func- 
tions are discussed. Various asymptotic approximations are derived 
and applied to analyze the basic structure of the two-dimensional 
Bessel functions as well as their nodal lines. 



1 Introduction 

Generalized Bessel functions depending on several variables have been intro- 
duced in 1915 for a finite [1] and also an infinite number of variables [2]. They 
have very similar properties as the ordinary Bessel functions but are much 
less familiar. More recently, however, they found an increasing number of 
applications in various areas of physics (see, e.g., [3-13]). The basic theory 
of generalized Bessel functions is described in a mongraph by Dattoli and 
Torre [14]. Our own interest into the properties of these functions is caused 
by our recent studies of quantum dynamics in periodic structures [11,12], in 
particular in studies of transport and dynamic localization [15]. 

In most cases these applications were restricted to the case of two vari- 
ables, u and V. Then the generalized Bessel functions J^''^{u,v) are labeled 
by three integer indices n,p, q. The special case (p, q) = (1, 2) has been con- 
sidered up to now almost exclusively [5,8,10,14,16-18]. Here we will analyze 
the two-dimensional Bessel functions J^''^{u,v) for general indices p and q 
(see [19] for a well written introduction to the case of infinite variables). 

We will derive the fundamental properties of the two-dimensional Bessel 
functions and analyze their basic structure for small and large arguments 

*Present address: Institut fiir Physik, Universitat Kassel, D-34109 Kassel, Germany 



1 



in the following sections. It will be seen that the two-dimensional Bessel 
functions show a rich oscillatory structure with regions of very different be- 
havior. We will analyze these structural features with special attention to 
the nodal lines which are of considerable importance for recent applications 
to localization phenomena in quantum dynamics [15]. 

2 Basic properties 

In this section we will collect the basic properties of the generalized Bessel 
functions J^''^{u,v) with integer indices n, p, q and two real arguments u, v. 
Most of the results in the literature (see, in particular, appendix B of [5] and 
chapter 2 of [14]) have been derived for the special case of J^'^('U, f) . 

2.1 Definition 

The two-dimensional Bessel functions can be defined by the generating func- 
tion 

oo 

^iiusinpt+vsingt) ^ ^ JP,1^U,v)e'"\ (1) 
n=— oo 

also known as a Jacobi- Anger expansion, or, somewhat more general, as 

oo 



ra=— oo 



Integration of ^ over t using J^J dt e"^* = 27r5„ o immediately leads to the 
integral representation 

Jl^\u,v) = — r dte'("^^'^P*+''^^'^^*-"*), (3) 
27r 

a generalization of the integral representation of the well known ordinary 
Bessel function 

JJx) = — r dte^("^'" (4) 
2vr 

From the properties of Fourier series we immediately find the bounds 

\Jl'%u,v)\<l and \JP'''{u,v)\ <1/V2 for n^O. (5) 
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As an immediate consequence of Q the integers p and q can be assumed 
to be coprime because J^'"^ vanishes otherwise or it can be reduced to such 
a coprime case. This is seen as follows (we assume that /x is an integer, 

f27r 

271- JMP>MCy t>) = / ^-i^qK^ sin fipt+v sin iiqt-nt) 







27r/i 1 A* i>2Trm i 

___ giCwsinps+ti sings— ns//i) \ ^ / ^if-u sinps+ti sin qs— ns//x) 

J27r(m-1) A' 

r 1 /'^''ds 

^ ^ i27r(m"l)n//j / ^if-u sinps+f sin gs— ns/^) 

Using 

r 

g-i27r(m-l)n/^ ^ I 

m=l 

and - for n/yU G Z - 



/i for n/yU G Z 
else 



(7) 



2vr J 1 

i{tismps+t) sm (js-ns/p) _ jP,q ( \ 

/X ~ /-/^^'''''^ 



one obtains 



4™^^^) = !.^""^"''^ forn^GZ 
[0 else. 

In the following, we will therefore assume that the integers p and q have no 
common divisor. 



2.2 Decomposition in terms of ordinary Bessel func- 
tions 

A representation in terms of ordinary Bessel functions can be derived from 
the integral representation 0. Inserting the generating function for the 
ordinary Bessel functions 

oo 

gixsin. _ J^^^^^ins ^^Q) 

n=—oo 
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for both s = pt and s = qt into 0, we obtain 

27r 

= ^J^(m) J,(t;)i- r"dte'[(^^'+"M. (11) 

The integral is only different from zero if n = fip + uq is satisfied. If p and q 
have no common divisor as assumed here, a solution (/i, z/) = (M, A^) of this 
Diophantine equation always exists and can be found systematically by, e.g., 
the Euclid algorithm [20]. Moreover there is an infinite number of solutions 
fj, = M — qk, u = N + pk, k = 0, ±1, ±2, ... because of 

n = pM + qN = pM + qN + pqk — pqk = p{M — qk) + q{N + pk) . 

We therefore have 

oo 

J^'\u,v) = Yl Jm-A^) JN+pkiv) , (12) 

fc=— oo 

where (M, A^) is an arbitrary solution of n = pM + qN. 
For the case p = 1 this reads (M = n and = 0) 

oo 

Jn''{%^)= J2 Jn-,k{u) Ju{v) . (13) 
fc=— oo 

In most of the previous applications one encounters the case q = 2 and in 
these cases one usually simplifies the notation by dropping the p, q indices, 
i.e. one defines 

Uu,v) = Jl'\u,v) . (14) 

2.3 Addition theorems 

The addition theorem 

oo 

Jr (^1 + «2, VI + V2)= Jn-k{u,, Vl) ^r(«2, V2) (15) 

fc=— oo 
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can be easily proved starting from Q using (jT)): 

JP.9(^y_j^ _|_ ^2 Vi + V2) = / fj^f-Qiiui sin pt+vi sin qt-nt) ^i{u2 sin pt+V2 sin qt) 

27r 

= — / ^tQKuisinpt+v,sinqt~nt) ^ jf'(M2, 
00 

= E •'r 

fc=— 00 
00 



V2 e 



1 

2^ 



(^U2,V2) / (^-i^qK"^! sin pt+vi sin qt-{n-k)t) 



k=—oo 



The Graf addition theorem for ordinary Bessel functions 

I _ n. 

X2 - Xi/t 



X2 - XiT 



Jn[9{Xl,X2;T)] 



(16) 



with 



g{xi,X2;T) = {xl + xj -XiX2(r + 1/r)) 



1/2 



(17) 



can also be generalized to the two-dimensional case, at least for p = 1 in the 
form 



+00 



t^jI'''{ui,vi)JI'1^{u2, 



V2 



e=-oo 



E 



£=-00 



U2 - Ui/t 


n — qi 
2 


V2 - Vi/t'^ 


_ U2 - UiT _ 




f 2 — ViTl 



■Jn~qe[9iui, U2; t)] Je[g{vi, V2; r'')] 



(see Dattoli et al. [14,16] for the special case q = 2). The generalized Graf 
addition theorem (fT^ can be derived in a straightforward calculation ex- 
pressing first the two-dimensional Bessel functions as a sum over ordinary 
ones (see eqn. (fT^ ) and using the Graf addition theorem (fT^ for ordinary 
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e=-oo 



Bessel functions: 

+00 

-00 

y^^r'^Je-qk{ui)JkiVl)Je+n-qj{u2)Jj{v2) 

Jk{Vl)Jj{v2) r^'' ^ T^' Ji>{Ui)Jn+q{k-j)+l'{u2) 



(19) 



3,k 



^ Jk{Vl)Jj{V2) r''^Jn+q{k-j){9{Ul, U2] t) ) 



E 
E 



U2 - Ui/t 
U2 - UiT 



n + q{k~j) 



U2 - Ui/t 
U2 - UiT 

U2 - Ui/t 

U2 - UiT 



71 — q£ 



Jn-qt{g{Ul, U2\ t)) ^ t''^ Jk{vi) Jl+k{v2) 



n — ql 



V2 - Vi/t'^ 
V2 — ViT'i 



Jn-qe{g{ui, U2] t)) J i,{g{vi, V2] t'')) 



with g{xi,X2;T) as defined in (fTTj) . 



2.4 Symmetries, special cases and numerical examples 

From the definition one verifies (by taking the complex conjugate and 
changing variables t — > —t) that the J^''^{u,v) are real valued and satisfy 

Jr(0,0) = 5„o- (20) 

The symmetry relations 

J'n'i^v) = Jr{v,u) (21) 

J^_!l{u,v) = JP''^{-u,-v) = JY{-v,-u) = J-^'-\u,v). (22) 

follow directly from the definition. For n = these equations imply the 
symmetries 



(23) 



A further direct result is a symmetry relation for even values of one of 
the p, g-indices, say q. Using Jn{—z) = (—1)" Jn{z) we get 

oo 

JP'1{-U,V)= ^ JM~qk{-u) JN+pk{v) 

k=—oo 
oo 

= (-1)^^ -hi-qkiu) Jm+M = i-ir J'n'%u, v) . (24) 



k= — oo 



The last equality holds because of (—1)" = (^—ly^^+i'^ = (—1)^'^ for q even 
and p odd. This symmetry implies 

JM(0, v) = JP'%0, v) =^ J^'^(0, v) = for n odd and q even. (25) 

If both upper indices are odd, their difference must be even. This leads to 
another symmetry relation 



J^'^(-U, -V) = ^ JM-qk{-u) JN+pk{-v) 
k=—oo 

oo 

= J2 {-lf'^''^^''-'^'JM-A^) Jn+Av) (26) 

fe=— oo 

= (-1)^^+^ Jl^\u, v) = i-ir Jr{u^ v) for p, q odd . 

Here the last equality is based on the fact that for odd p, g-indices, p = 2j + 1 
and g = 2/c + 1, we have n = pM + qN = M + N + 2{j + k). 

In the case p = q the two-dimensional Bessel functions simplify and reduce 
to ordinary Bessel functions if n is an integer multiple of p: 

1 r+n 1 p+pn j 
JP'Pfll y) = / (^fQ^iiu+v) sm pt-nt) _ / ^i{(u+v) sm s-ns/p) 

Jn/p{u + v) for n/p e N , . 

else. 

Another relation between the generalized and ordinary Bessel functions can 
be observed if the index n is a multiple of one of the upper indices, e.g. 
n = mq, m integer. Then we get 

JP^I{{],V) = — f^tQiivsinqt-mgt) ^ / ^^v sin s-ms) 

2^ J-n 27rq' J_qj^ 

= -1 r"dse'(''^''^™) = J„(i;) (28) 
27r 

and, as a special case, 

JP'\0,v) = Uv). (29) 

According to the two-dimensional Bessel function Jn'^{0, v) is reduced 
to an ordinary Bessel function Jm{v) for n = mq. Otherwise the function 
vanishes on the f-axis as can easily be seen from ()12j) : 

oo oo 
J^'^(0,t;)= ^ JM-qk{0) JN+pk{v)= ^ SM,qkJN+pk{v) = , (30) 
k=—oo k=—oo 




-20-10 10 20 
u 



-20-10 10 20 
u 



20 
10 

> 



-10 
-20 



-20-10 10 20 
u 



Figure 1: Color map of the two-dimensional Bessel function J^''^{u,v) for 
n = 0,n = l,n = 2 (from left to right). 
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Figure 2: Color map of the two-dimensional Bessel function Jl^^{u,v) for 
n = 0, n = l, n = 2 (from left to right). 



if M is not a multiple of q or, equivalently, n = pM + qN is not an integer 
multiple of q. We therefore have 



JP''?(o, v) = a n ^ mq, m e Z 



(31) 



generalization of ()25j) . 

Let us look at a few examples of two-dimensional Bessel functions calcu- 
lated numerically using the representation (fT^ in terms of ordinary Bessel 
functions (similar graphs can be found, e.g., in [14,17]). Figures [T] and |21 
show color maps of J^'^(m, v) and Jn^{u, v) for n = 0, n = 1 and n = 2 using 
a re-normalization to unit maximum in each case (the regions of positive 
values are colored red, of negative ones blue). 

For J^'^ in figure ^ we see that the symmetry relations 



J^'\-u,~v) = 4'\u,v) , J'^'\-u,v) = i-ir 4'\u,v) (32) 

are satisfied (compare eqs. and ()24j) ). The last relation implies a nodal 
line for = 1 along the f-axis, jI''^{0,v) = 0. For the case J^'^{—u,—v) 
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in figureElwe have the symmetry J^'^{—u, —v) = (—1)" Jn'^{u,v) (compare 
eq. m)- 

2.5 Sum rules and Kapteyn series 

The simple sum rule 

oo 

J2j'n'M = l (33) 

?i=— oo 

is a direct consequence of the generating function for t = 0. A variety of 
sum rules for special cases can be obtained by choosing t in ((T)) appropriately. 
E.g. for the important special case {p,q) = (1,2) another sum rule is found 
by setting t = n/2: 

oo 

^ i-Jl'\u,v) = e-. (34) 

n=— oo 

Similar sum rules can be obtained for other special cases. Another sum rule, 

oo 

J2 {jr{u,v)Y = i (35) 

fc=— oo 

follows from the addition theorem (fTH|l for the special case n = 0, ui = —U2 = 
u, vi = —V2 = V using (1201) • 

We furthermore note without proof the Kapteyn type series [21,22] 

oo ^ 

7 JI!''^(nu,nv) = , b-ul + |gt>| < 1 . (36) 

^-^ 1 — pu — qv 

n=— oo 

2.6 Further generalizations 

As already stated in the introduction, the number of variables in the Bessel 
function can be extended. Different types of generalizations are, however, 
also possible. Modified higher dimensional Bessel functions can be con- 
structed, e.g. by replacing one of the ordinary Bessel functions in ()12|1 by 
a modified one [9]. In addition, two- variable, one-parameter Bessel func- 
tions [9,23] can be defined generalization of ()12j) : 

oo 
fc=— oo 
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Figure 3: Color map of the two-dimensional, one-parameter Bessel function 
J^'^(u, f ; i) for n = 0,n = l,n = 2 (from left to right). The figures show the 
real part. 



(Let us recall that {M,N) are arbitrary solutions of n = pM + qN.) Here 
again we find J^'''(0, 0; r) = 6no ■ 

In particular the case r = e"^ is of interest [9, 15] for applications in 
physics. We confine ourselves to the most important case p = 1, i.e. 

oo 

Jl'\u,v-e'') = Yl Jn-A^) Jk{v)e''' . (38) 

fc=— oo 

Following the lines in the derivations above, one can easily show that these 
functions are generated by 



oo 

^i{u sin t+v s'm{qt+5)) 



J^'^,v;^')^\ (39) 



which leads to the integral representation 

J^'''{u,v;e'^) = — r dte^("''''*+^'"^(^*+'^)-"*) . (40) 
27r J^^ 

The generalized Bessel functions J^'^^u, v; e^^) satisfy most of the properties 
of the J^''^{u,v) , as for example the bounds (0), the addition theorem (fT3|) 
and the sum rules ()33|1 and (j35|) . These function are, however, complex 
valued. 

Figure El shows the real part of the Bessel functions J^'^(m, f;i) for n = 
0, 1, 2. A comparison with figure [T] shows that the structure of the functions 
is strongly altered by the angle parameter S. 
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2.7 Differential equations and recurrence relations 

Finally, we will derive recurrence relations for the J'fl''^{u, v) and their deriva- 
tives. Differentiating the generating function (j21) with respect to u leads 
to 

i(.--.-)et(-^--"^) + i(^^--"^) 

oo oo 

= i '^n'i^^, V) (^"^" - ^"-") = E ^uJ'n'i^, V) . (41) 

n=— oo n=— oo 

Equating the coefficients of z", we find 

2 d^JP''^iu, v) = Jl\{u, v) - J:Uu, v) (42) 

and similarly 

29, (w, v) = J^^,(w, v) - J^^^iu, v) . (43) 

If we differentiate (j21) with respect to z and compare the coefficients, we find 
the recurrence equation 

pu (J^^,(w, v) + J::'Uu, v)) +qv (J^\(w, v) + J^'^u, v))=2nJP^'^{u, v) . (44) 

These are generalizations of the relations derived by Reiss [5] for the case 

Similarly one can show that the derivatives of the generalized Bessel func- 
tions (PH|1 are given by 

2dX''i^^, v, e'') = e^'jl'l^iu, v; e'') - e-'V^^,(w, v; e''). (45) 

Using these relations one can show that the two-dimensional Bessel func- 
tions solve a variety of linear partial differential equations, depending on 
their indices (p, g). These differential equation can be constructed system- 
atically by adding up derivatives of different order such that the different 
terms J^'*^ cancel each other for every value of n. For example the ordinary 
two-dimensional Bessel functions with q = p solve the wave equation 

{dl-d',)j'/iu,v) = 0, (46) 

while the generalized Bessel functions with {p,q) = (1,2) and 6 = n/2 solve 
the time-dependent Schrodinger equation [17, 24] 

id.j'Jiu, V- i) = {-2dl - 1) Jl^\u, V- i). (47) 
11 



Furthermore, a repeated application of the differentiation rules fH^ . (jlHj) 
and the recurrence relations (jl^ leads to the coupled differential equations 

[{pudu + qvd^f + {p^udu + q^vdy) + p^n^ + q^v"^ - n?] Jl;'^{u, v) 
= -pquv ( J^V.(^' ^) + Jnlr>-,i^. v) - 2 (w, v)) (48) 

for arbitrary indices {p,q)- Except from the right-hand side this equation 
is structurally similar to the defining differential equation of the ordinary 
one-dimensional Bessel functions. Equations (jlHj) can be decoupled for q = 
up, z/ G Z by applying again (jUJ. For (p, q) = (1, ±1) this yields 

[{udu + vd^f + [udu + vd,) + {u± vf - n^] Jl^^\u, v) = 0, (49) 

while the respective calculations for other values of q lead to more compli- 
cated results. 



3 Polynomial expansion for small arguments 

We will first analyze the regime of small arguments u and v and derive a 
leading order poljTiomial expansion. Following Wasiljeff [4], we expand the 
u, f-dependent part of the exponential function in (jHl) in a Taylor series: 

— TV 

, oo ^ 

27r^2^fc! J ^ ' ^ ' 

Using the polynomial formula 

{a + h + c + dy =3\Y, , .TTT^TTTT, (52) 
^ — ^a,/3,o-,c a\ p\ a\ Q\ 

where the primed sum runs over all indices with j = a + (3 + a + (, one 
obtains after rearranging terms and carrying out the integration the series 
expansion 



oo 



1 " „,a+/3„.CT+C 

Here the double-primed sum includes all nonnegative integers with 

j = a + f3 + a + C and n = {a - f3)p + (a - Qq . (54) 
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The sum can be transformed into a more convenient form by introducing 
i = a + (3, 2f = i + a — l3,m = a + ( and 2g = m + a — Q. After some 
elementary algebra, this yields 

^rK-)= E -tr'%i (55) 

with ^ 

m-fmrn-g)\- ^^^^ 

g=0,...,+m 
n=p{2f-l)+q{2g~m) 

The lowest order approximation in this expansion can be found in explicit 
form for the case p = 1. Then the lowest order term in is given by^ 

(n-g[|J, 0, [^J, 0) if I ^ N and 2gmod(n, g) <q + l 

(«, (5, a,C) = { (0, -n + g [H] , pn] , q) if f ^ N and 2g mod(n, g) > g + 1 

(0,0,^,0) if^GN. (57) 

In most cases it is given by a single term, as for example in 

22 ' -4^-5-/ 22 

1 2 

Note that, because n = 3 is not a multiple of g = 2, the Bessel function Jg' 
is identically equal to zero on the f-axis (see eq. (jSH)), however Jg'^ and J^'^ 
do not vanish on the u-axis, where we have J^'^(m,0) = Jniu) ~ {u/2Y/n\ 
(see eq. p9|l V 

In certain cases more terms of the same minimum order j appear. This 
happens if ^ ^ N and 2gmod(n, g) = g + 1. One can easily check that this 
requires that g is odd, g = 2z/ + 1, and n = nq + u + 1 with z/, G N. In this 
case, the lowest order approximation reads 

This yields the (approximate) nodal line 

V = u (d(J) 

for small u and v. As an example, we note g = 5 and n = 23, i.e. u = 2 and 
fi = 4 and therefore 

15, , u'^v^ /M f 



J3'^(u, t;) ~ — Mt; , J]'^(m, w) ~ — t;^ (5^ 



J2T(n, v) ~ — — - + - . (61) 
23 V . ; 272!4! V3 5/ ^ ' 



"'^Notation: [zj is the largest integer < x, \x] is the smahest integer > x and 
mod(n,g) = f -LfJ. 
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4 Asymptotic approximations 

In the examples of two-dimensional Bessel functions J^''^{u,v) shown in fig- 
ures [T] and El for {p,q) = (1, 2) and {p,q) = (1,3), respectively, one observes a 
rich oscillatory structure which will be analyzed in the following. The skele- 
ton of this structure and valuable approximations can be obtained asymp- 
totically by means of the stationary phase approximation 



The sum extends over all contributing real stationary points and the ± sign is 
chosen so that ±g"(ts) is positive (see, e.g., [25] for more details). Previous 
studies of asymptotic approximations for two-dimensional Bessel functions 
[3, 10,26] have been restricted to the case p = 1, q = 2 and special regions of 
the index n and arguments u, v. 

Information about the oscillatory structure of the multivariable Bessel 
functions J^''^{u, v) can be obtained from asymptotic approximations for large 
arguments and/or large indices. We will consider three of the large number 
of possible limits: the case when both arguments u and v are large, whereas 
n remains fixed, and the case where one argument, v, and the index n are 
large for fixed value of the argument u. Finally we will consider the limit 
where both variables as well as the index n are large. In all cases we assume 
small fixed values of the indices p and q. 

4.1 Basic structure for large arguments u and v 

We base our analysis on the integral representation 



jP'i(u,v) = — dte'("^(*)""*) , (f){t) =u sin pt + v sin qt. (63) 



Here we will consider the asymptotic limit of large arguments u and v as- 
suming that n is fixed, i.e. we identify g{t) = (j){t) in an application of (jU^ . 
The condition 



determines the stationary points ts (note that there are always pairs of such 
stationary points with different sign due to the symmetry of the cosine- 
function). The integral (jU^ is then approximately given by 



TT 




(62) 




(p'if) = pu cos pt + qv cos qt = 



(64) 




(65) 
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where the ±-sign is given by the sign of (j)"{tj). The main contribution to 
the sum is provided by real-valued stationary points, complex points lead 
to exponentially decaying terms. At the points where two stationary points 
coalesce when the arguments u and v are varied, the second derivative 

(t)"{t) = —p'^usinpt — q^vsinqt (66) 

vanishes and the approximation diverges. Crossing these bifurcation points, 
the function changes its character. In the present case, the bifurcations are 
determined by the simultaneous solution of 

pu cos ptg = —qv cos qts and p'^u sinptg = —q'^v sin qts ■ (67) 

This can be most easily satisfied if both sides of one of the two equations are 
equal to zero. We distinguish two cases: 

Case (i) sinptg = — smqtg = : For coprime p and q, this implies = or 
tg = IT and therefore (from pucosptg = —qv cos qtg) we have pu = —qv or 
pu = —{—ly^'^qv, respectively. We therefore obtain the bifurcation lines 

V = ±pu/q if one of the p, q is even (68) 

V = —pu/q else . (69) 

Case (ii) cos ptg = cos qtg = : This implies pt^ = vr/2+j7r and qtg = 7r/2+A;7r 
with integer j and k or {2k + l)p = {2j + l)q and (for coprime p and q) 
p = 2j + 1 and q = 2k + 1 . With sin ptg = sin(7r/2 + jn) = (—1)"' and 
sinqts = sin(7r/2 + kn) = (— l)'^ the second condition in (j67j) leads to 

V = -{-iy+''up'^/q'^ p = 2j + l,q = 2k + l. (70) 

The examples in figure[T]show Bessel functions Jn'^iu, v) for various values 
of n. Here q is even and from equation ()69|1 we find the bifurcation lines 

v = ±u/2. (71) 

We observe that the structure of the Bessel functions changes if one crosses 
these lines. In the left and right sectors, we have only two stationary points, 
±ti, whereas in the upper and lower sectors we have four stationary points, 
±ti and ±)f:2, and consequently a richer interference pattern. This will be 
analyzed in more detail below. 

For the two-dimensional Bessel function J^'^(m, w) displayed in fig. Elboth 
upper indices are odd and the bifurcation lines are given by eqs. ()69p and 
dZOl): 

f = — u/3 and v = u/9 . (72) 
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The qualitative difference to the behavior of J^''^{u,v) in fig. ^ is obvious. 

Let us now analyze the function J^'^(u,f) in more detail working out 
explicitly the stationary phase approximation. In view of the symmetry 
Jl;'^{u,v) = j]_'^{—u,—v) (j22I) we can assume > in the following for 
simplicity. The stationary phase condition 

ucost = -2v cos2t = -2v (2cos^t - 1) (73) 

can be solved for c = cost with solution 

c± = ^{-u/v± Viu/vy + 32 ) . (74) 

In the region —2 < u/v < +oo the necessary condition |c±| < 1 is met by 
c+ and vice versa by c_ in the region — oo < u/v < +2. Note that in the 
interval —2 < u/v < +2 both solutions fulfill |c-i-| < 1. With t± = arccos c± 
and sint-i- = ±a/1 — we arrive at 



u smt± + V sin 2t± = ±{u + 2vc±) \/l — c| (75) 



and with the definitions 



" - -Msint± - 4t;sin2t± = =f(m + 8fc±) Wl - 4 (76) 



, ^ , -rTTTT cos ((/)+ — n arccos c+ — 7) for — 2t> < m , , 
F+{u,v) = <{ V ^ (77) 

else 



cos 



+ n arccos c_ — f ) for m < +2v 



F.iu,v) = <; V^l^-I ' 4^ (7J 

else 



the final result can be written as 

Jl'\u,v)c:^F+{u,v) + F.{u,v). (79) 

Here one should be aware of the fact that in the region |n| ^ 2\v\ both of 
the terms F±{u,v) provide a non- vanishing contribution. 

This so-called 'primitive' stationary phase approximation diverges at the 
bifurcation lines v = ±2u. If desired, it can be improved by taking complex 
stationary points into account and by taming the divergences by uniformiza- 
tion methods. 

In the limit |m| — 00 the asymptotic approximation ()79p simplifies dras- 
tically. Only contributes for u > (F„ for m < 0) and with t± = ±7r/2, 
0± = 0± = we find 



Jn'\u, v)c^J—- COS {u-n---), (80) 
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Figure 4: Two-dimensional Bessel functions J^''^{u, v) with f = 10 for n = 
(left) and n = 1 (right) as a function of u (full curve) in comparison with the 
stationary phase approximation ()79|) (open circles). 



which agrees with the well-known asymptotic approximation of the ordinary 
Bessel function J„(m) for large arguments [27]. In the alternative limit v oo 
both terms, F+ and F_ contribute. With c± = —-^ ± (j)± = v ± -j^ and 
(f)'^ = ±4f we obtain 

/ 9 f + cos (v — (n + 1)^) cos n even 
J^'^(.,.)^ . ^ ^^^^^ . f (81) 




— sin (t> — (n + 1)^) sin odd , 



as already derived in [14]. 

As an illustration of the asymptotic formula ()79p . fig. ^ shows the two- 
dimensional Bessel function J^''^{u,v) in comparison with the stationary 
phase approximation for v = 10 and n = 0, 1. With the exception of 
the vicinity of the divergences at u = ±2v = ±20, the agreement is ex- 
cellent. This approximation can be used in oder to determine the nodal lines 
of the two-dimensional Bessel functions which is of interest for applications 
in physics [15]. 



4.2 Large argument v and large index n 

In this section we will consider the regime where the index n and one of the 
arguments, e.g. v, are large. In view of ()22|) . we can assume n > 0. Following 
the analysis applied to the special case J^'^(M,f) by Reiss and Krainov [10], 
we separate the integral representation (jHJ into a fast, e^^*^*-*, and a slowly 
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oscillating part: 

TT 

JP^i(u, v) = ^ [ dt e^'^^i'^PV^W , g{t) = v sin qt - nt , (82) 
27r J 

— TT 

and evaluate the integral approximately by the method of stationary phase 
or the saddle point method if the stationary points are complex valued (for 
details see, e.g., [25]). 

The stationary points tg of g{t) are obtained from g'{ts) = as 

n 

cos qtg — — (83) 
qv 

with real valued solutions for n < q\v\ and complex solutions otherwise. We 
discuss these two cases separately. 

(i) For n < q\v\ the stationary points are 

1 n 

tf = ±(to + 27rs/q) , s = 0, 1, 2, . . . , to ^ - arccos — , (84) 

q qv 

i.e. a finite number in the interval — tt < tf < tt. 
With sin qtf — ± sin qtg we have 

g{tf) = vsinqtf — ntf = ±vsmto^n{to + 2ns/q)] (85) 

9"i^f) — —q^v sin qtg = ^q^v sin qto . (86) 

and, using sinptf — ±sinp(io + 27rs/g), the final result is 

JP'1(U V) ~ / 1, ±i[«smp(to+^)+i'smgto-n(to + ^)-f ] 

n \ 1 J V 2irq'^vsinqto ^ 

- S cos {u sinp(to + ^) + ^ sin qto - n{to + ^) - f ] • (87) 

We will work out the case p — 1 and g = 2 in more detail. Here we find four 
stationary points ±(io + stt) with s — and —1 and therefore 

J^'^u, v) ~ ^ 2nvL2to { COS [ M sin to + v sin 2to - nto - ^] 

+(-1)" cos [ - u sin to + v sin 2to - "-^o - ^ ] } (88) 
+ cos(-u sin to) cos(f sin 2to — nto ~ j) n even 



TTU sin 2to 



■ sin(ti sin to) sin{v sin 2to — nto — f ) n odd 
18 



with 



1 ?1 / Ti^ /I Ti 

to = -arccos— , sin2to = yi-^, sinto = ^--—. (89) 

In comparison with the semiclassical approximation derived in section 14.11 
the result ()88p agrees approximately with (f?^ also for small values of n, as 
for example Jq''^{u,v) shown in figure 0] for v = 10. Equation (jHHj) misses 
however the structural transition at \u\ = 2v and cannot describe the region 
\u\ > 2v. 

(ii) For n > q\v\ the stationary points ()84|) are complex, tg = Xg + iy with 
real part 

X = ( ^^""/^ ' > ° s = ±1 ±2 f90) 

\{2s + l)n/q ,v<0 U, ±i, ±2, . . . [JO) 

with — TT < Xg < +TT- The imaginary part is the same for all s: 

y± = ±-arccosh( — — - ) . (91) 
q \q\v\J 

The integral is approximately carried out by the saddle point integration, 
where the integration path is deformed to a steepest decent curve passing 
through the saddle points [25]: 



TT 



dt hit) e^^W ^ Yl V ^ M^s) e'^(*^) . (92) 

s 

The second derivative is 

ig"{ts) = —ivq'^ sin qts = —q^vsinhqy (93) 

and the conditions for the integration path [25] can only be satisfied for the 
saddle points in the upper (lower) complex plane for f > [v < 0). With 

ig{ts) = i{v sin qtg — ntg) = —v sinh qy — inxg — ny (94) 

we obtain the result 

-|vsinhijj/|-n|j/| 

jp^i[u, v) = V e-inx^ei^^i'^pi^^+iy)^ (95) 

y 27rg2|f sinh gy| 

Let us again consider the special case p = 1, g = 2 in more detail. Two 
saddle points contribute [xg = ±7r/2 with ?/_ for w < or = and tt with 
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y+ for t> > 0) and equation simplifies. 
For V < we find 

JI''^{u,v) = , COS (m coshu —n — ] (96) 

with 



1 n hn? n 1 

V = - arccosh— — , sinn2'u = \ -— — 1 , coshy = a / — — - H — (97) 
^ 2 2|v| ' ^ V 4f 2 ' ^ \j 4|f| 2 ^ ^ 

in agreement with the resuh derived in [10]. 

For V > the resulting approximation is non-oscillatory: 

e"''*'''''^^^""^ f cosh(Msinhv) n even 
-T':^^^' = — ^ ^ ^ (98) 

V47rV2t; sinh 2?/ [ sinh(Msinht/) n odd 

with 



1 n n 1 
y = - arccosh — , sinhu = \ . (99) 



Note that both asymptotic approximations (|88|) and (j9H|) satisfy the symme- 
try relation Jn{-u,v) = (-1)" Jn(M.w) (cf. eq. ^ ). 

Figure El demonstrates the quality of the asymptotic approximation for 
n = 30 and v = 64 (case (i)) or f = — 12 (case (ii)). Reasonable agreement 
is observed for |m| < 10. These simple approximations get worse in the 
vicinity of n = q\v\ where they diverge. A finite result can be obtained 
using an appropriate uniformization technique, in the present case an Bessel 
uniformization, e.g. a mapping onto an (ordinary) Bessel function [28] (see 
also [26] for an alternative method). 

The asymptotic approximations (jHH)) and provide explicit estimates 
for the nodal lines of J^'^('U, f). For n < 2\v\ we find 



/l ^ I (2j + l)f n even 
and for n > 2 If I we have for f < zeros at 



1 Ti TT 

uM-- — ={2j + n)- , J = 0, ±1, ±2,... . (101) 
These results are, of course, in agreement with the zeros observed in fig. El 
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Figure 5: Comparison of the asymptotic approximation (j96|l (open circles) 
with the exact two-dimensional Bessel function J^''^{u,v) (full line) for n = 
30. Left: Case (i) for v = QA using eq. (jHHl). Middle: Case (ii) with v = -12 
using eq. (|^. Right: Case (ii) with v = +12 using eq. 



4.3 Large arguments u, v and large index n 

As an example of the structure of the two-dimensional Bessel functions for 
large indices, fig. El shows Jl'^{u,v) for n = 29 and n = 30. These functions 
look quite similar, they are clearly distinguished, however, by their symmetry 
property J^'^(— M,f) = (— 1)" J^'^(m, f ) (see eq. (j211)), i-e. J^l^ is even and 
is odd with respect to a reflection u —>■ —u. Therefore vanishes 
on the f-axis, ^29^(0,^) = 0. The function Jgg^ is symmetric on the f-axis: 
J^'^(0, — f ) = J^'^(0,f) (see eq. (jSI}), despite of the apparent asymmetry 
with respect to the reflection v —v. 

In additions to the oscillatory pattern in the four sectors, we observe a 
region close to the center where the values of the Bessel functions are small. 
This pattern can again be explained by a consideration of the asymptotic 
limit where both arguments and the index n are large using 

g{t) = usinpt + vsmqt — nt (102) 

in the stationary phase approximation (jH^ . The stationary points ts are 
determined by 

9 {is) = W cospts + qv cos qtg — n = . (103) 

The zeros of the second derivative 

g"{t) = —p^usinpt — q^vsinqt (104) 

appearing in the denominator of ()62|) determine the bifurcation set of these 
solutions. 
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Figure 6: Color map of the two-dimensional Bessel functions J^''^{u,v) for 
n = 29 (left) and n = 30 (right). Note the different symmetries of these func- 
tions. Their overall structure can be explained by means of the bifurcation 
set shown in fig. d 



Restricting ourselves again to the case {p,q) = (1,2) these equations 
simplify and can be solved in closed form: 



with solutions 



C± = costs 



u COS ts + 2v COS 2ts — n = 



u 



± 



u 



2 1 n 

+ o + ^ 
2 Av 



(105) 



(106) 



(here again each solution c± implies two stationary points ts because of the 
symmetry of the cosine function). The bifurcation set - the skeleton of the 
Bessel function - is found when 

g"{ts) = -u sin ts - Av sin 2t, = (107) 

is satisfied in addition to ()105|) . Eliminating ts we find the solutions 

v = {n±u)/2 (108) 

(note that for \u\ ^ n these straight lines agree with the ones stated above 
in eq. (f7T|) ) and the ellipse 



lQ{v + n/Af 



2v? 



1 



(109) 
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Figure 7: Bifurcation curves of the stationary points for the two-dimensional 
Bessel function jIq{u,v). 

centered at (m, f) = (0,— n/4) with half axes \/2n and n/4. Inside this 
ellipse the stationary points ()106p are complex, outside they are real. A brief 
calculation furthermore shows that the straight lines ()108|1 are tangential to 
the ellipse This bifurcation set is shown in fig. [3 In the upper sector 

(I) between the bifurcation lines we have —1 < c± < +1, as well as in the 
lower sector (II) outside the ellipse. Hence we have four real solutions ts in 
these regions and a corresponding oscillatory pattern. In the right sector 

(IV) two of these solutions become complex because of |c+| > 1 and simi- 
larly in the left hand sector (III) with |c_| > 1. In the triangular segment 

(V) in the lower sector above the ellipse, we have |c+| > 1 and |c_| > 1 
and therefore no real stationary points. In the elliptic region with complex 
valued stationary points, the Bessel function is damped but still oscillatory. 
An example is shown in fig. El (right hand side) which shows a cut through 
the Bessel function jIq{u,v) shown in fig. IHlfor v = —12 close to the elliptic 
bifurcation curve. A cut at f = 64 (left hand side) shows the oscillations 
in region (I). Note that the semiclassical approximations shown in fig. El are 
the simplified versions developed in section 14.21 A more refined semiclas- 
sical treatment along the lines discussed above will provide a much better 
agreement for larger values of u (compare also the treatment in [5]). 
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